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In this paper we study an archetypical scenario in which repulsive Casimir-Polder forces between 
an atom or molecule and two macroscopic bodies can be achieved. This is an extension of previous 
studies of the interaction between a polarizable atom and a wedge, in which repulsion occurs if the 
atom is sufficiently anisotropic and close enough to the symmetry plane of the wedge. A similar 
repulsion occurs if such an atom passes a thin cylinder or a wire. An obvious extension is to compute 
the interaction between such an atom and two facing wedges, which includes as a special case the 
interaction of an atom with a conducting screen possessing a slit, or between two parallel wires. 

To this end we further extend the electromagnetic multiple-scattering formalism for three-body 
interactions. To test this machinery we reinvestigate the interaction of a polarizable atom between 
two parallel conducting plates. In that case, three-body effects are shown to be small, and are 
dominated by three- and four-scattering terms. The atom-wedge calculation is illustrated by an 
analogous scalar situation, described in the Appendix. The wedge-wedge-atom geometry is difficult 
to analyze because this is a scale-free problem. But it is not so hard to investigate the three-body 
corrections to the interaction between an anisotropic atom or nanoparticle and a pair of parallel 
conducting cylinders, and show that the three-body effects are very small and do not affect the 
Casimir-Polder repulsion at large distances between the cylinders. Finally, we consider whether 
such highly anisotropic atoms needed for repulsion are practically realizable. Since this appears 
rather difficult to accomplish, it may be more feasible to observe such effects with highly anisotropic 
nanoparticles. 

PACS numbers: 42.50.Lc, 32.10.Dk, 12.20.-m, 03.70.+k 


I. INTRODUCTION 

It is quite remarkable that after nearly seven decades, 
interest in the so-called Casimir effect l] remains so high. 
There have been many theoretical and experimental de¬ 
velopments in the last few years. For a review of the 
status of quantum vacuum energy phenomena in general, 
the reader is referred to the volume edited by Dalvit et 
al. 0. 

One of the hottest topics in the field is the subject of 
repulsive Casimir effects. This could have a major im¬ 
pact in nanotechnology, where at distances well below 


•Electronic address: milton@nhn.ou.edu 
1 Electronic address: ekabalo@gmail.com 
1 Electronic address: prachi@nhn.ou.edu 
§ Electronic address: nima.pourtolami@gmail.com 
^Electronic address: iver.h.brevik@ntnu.no 
•‘Electronic address: simen.a.ellingsen@ntnu.no 
^Electronic address: stefan.buhmann@physik.uni-freiburg.de 
^Electronic address: stefan.scheel@uni-rostock.de 


1 /im, Casimir forces can play a dominating role. Re¬ 
pulsion can occur between electric and magnetic bodies, 
between electric bodies separated by a medium with an 
intermediate value of permittivity, or purely due to the 
geometry of the two bodies. The geometric repulsion 
that was demonstrated in Ref. [s] was based on numer¬ 
ical methods. A great deal of effort has been given to 
understanding the underlying analytical structure of re¬ 
pulsive effects 0-01 • A brief review, with references, is 
provided in Ref. [7]. For further work on repulsion see 
Refs. dm. 

In this paper, we will further develop the multiple¬ 
scattering approach to include three-body effects, which 
was introduced in the scalar context in Refs. [1 ll4l3l |. 
in particular in the context of Casimir repulsion. The 
electromagnetic formulation is in turn based on Green’s 
dyadics, which have a long history. The Green’s dyadic 
approach to computing the Casimir effect was first pro¬ 
posed in Refs, [lj, [l]|. This was a tensorial generaliza¬ 
tion of the scalar Green’s function variational approach 
Schwinger had given a few years earlier [l6|. All of this 
was in the direct line of evolution to what is now referred 
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to as the multiple scattering method. 


II. THREE-BODY CASIMIR ENERGY 


Although it is well-appreciated that Casimir forces 
are not additive, most work on such interactions has 
concentrated on forces between two bodies. But the 
multiple-scattering formalism is easily generalized to in¬ 
clude three-body interactions. In Sec. Hof this paper we 
develop the three-body formalism for the electromagnetic 
case. This is relevant to computing the force between a 
polarizable atom and two co-planar half planes, forming 
a slit, for example. As a simple illustration of the formal¬ 
ism, we re-examine, in Sec. mu the interaction between 
an atom and two parallel conducting plates, first consid¬ 
ered by Barton [13] • Naturally, the two-body forces dom¬ 
inate when the atom is near one or the other of the plates, 
and three-body effects become significant at the several 
percent level only when the atom is roughly equidistant 
from the two plates, but then the force is quite small. 
In Ref. [IS;] we examined the nonmonotonicity that can 
arise when two polarizable atoms are near each other 
and close to a conducting plate. These are reminiscent 
of effects seen between two macroscopic objects and a 
wall 0,0; our work generalized that given in Ref. [0 ■ 
We broke up the three-body terms into three- and four- 
scattering contributions; although both are comparable 
at short distances, as expected, the former dominate for 
atoms far from the plate. In order to work out the three- 
body effects for an atom interacting with two half-planes, 
constituting a slit in a conducting plane, or more gener¬ 
ally, facing wedges, in Sec. IIV Al we work out the scat¬ 
tering matrix for a single wedge, and then, in Sec. IIVB1 
apply that to recalculate the two-body repulsion found 
in Ref. 3. Because of the complexity of the calculation, 
a scalar analog is also considered in the Appendix. We 
then go on to consider three-body effects between a po¬ 
larizable atom and a pair of parallel conducting cylinders 
in Sec. El In Ref. @ we showed that for an anisotropic 
atom moving along a line perpendicular to but not in¬ 
tersecting a perfectly conducting cylinder, and polariz¬ 
able along that same line, a repulsive force occurs near 
the cylinder, provided the distance of closest approach is 
sufficiently large compared to the radius of the cylinder. 
(The same does not occur for a sphere.) Here we consider 
the three-body effects due to a second cylinder parallel 
to the first, so the pair forms an aperture, perpendicu¬ 
lar to which the anisotropic atom moves. We adapt our 
formalism to this case, where a multipole expansion is 
also possible, and show that when the distance between 
the cylinders is sufficiently large, repulsion is not affected 
by the three-body corrections, since the latter are very 
small. Finally, in Sec. EH we present a calculation that 
suggests that highly anisotropic atoms, necessary to ex¬ 
hibit the repulsive effects we are considering, may be be¬ 
yond reach. Therefore, as in the numerical calculations 
of Ref. Q, it may be more appropriate to consider the 
interaction with highly anisotropic conducting nanopar¬ 
ticles, such as needles, as suggested in related work on 
negative Casimir-Polder entropy [22}. 


The multiple scattering formulation has proved excep¬ 
tionally useful in computing Casimir energies for com¬ 
plex configurations. It is usually presented in terms of 
potentials, where the potential stands in for the devia¬ 
tion of the permittivity from its vacuum value, for in¬ 
stance. Here, however, we wish from the outset to con¬ 
sider perfect conductors, so we give the formulation en¬ 
tirely in terms of scattering matrices. In particular, we 
wish to analyze three-body effects. The formalism we 
apply here appears in many places; recent examples are 
Refs. 0,0. We use natural units, with H = c = 1, and 
Heaviside-Lorentz electromagnetic units, except for the 
definition of the polarizability. 

The quantum vacuum energy, with the bulk vacuum 
energy subtracted, is in general given by 

E = ^TrlnTTo 1 , (2.1) 

where the Tr symbol represents a trace over tensor indices 
as well as spatial coordinates. Here To is the free Green’s 
dyadic, which, for a given frequency w, can be written as 

T 0 (r, r') = (1 uj 2 + VV)G 0 (|r - r'), (2.2) 

in terms of the free Helmholtz Green’s function 

„i\ui\R 

Go(R) = R = |r —r'|. (2.3) 

The full Green’s dyadic T satisfies the same differential 
equation as the free Green’s dyadic, 

r-'r = i, (2.4) 

where 

T" 1 = -^VxVx-l = [VV-(V 2 +u 2 )l] . (2.5) 

Here we have adopted a matrix notation for both the 
tensor indices and the spatial coordinates, so 

1 = l(5(r — r'), (2.6) 

where on the right 1 refers to the tensor indices only. The 
conducting surfaces S appear through boundary condi¬ 
tions on the Green’s dyadic, 


n x r 


= 0 , 


(2.7) 


where n is the outward normal to the surface at the point 
in question. We may define the scattering matrix T by 


so that 


r = r n 


T = -To 1 


r 0 Tr 


0 , 


r^riv 1 . 


( 2 . 8 ) 

(2.9) 
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Now we turn to the quantum interaction of three bod¬ 
ies. It seems easiest to start with the situation where the 
bodies may be described by potentials Vj, i = 1, 2, 3, 
and then write the result in a form in which only the T 
operators appear, so it applies to the conducting bound¬ 
ary problem, defined by Eq. m- The total potential is 
V = Vi - 1 - V 2 TV 3 , and the vacuum energy is given by 
the trace-log of 

rrv 1 = (i-r 0 v)-\ (2.10) 


or 

E = -^Trln(l - r 0 V). ( 2 . 11 ) 

Now it is easy to see that 

1 - iyvi + v 2 + v 3 ) = (1 - t 0 v 1 - r 0 v 2 ) 
x [1 - (1 - r 1 v 2 )- 1 r 1 v 2 r 1 v 3 (i - r^y 1 ] 
x(i-r 0 v 1 )- 1 (i-r 0 v 1 -r 0 v 3 ). (2.12) 

Here we have introduced the Green’s dyadic belonging to 
potential i alone, 

r* = (1 - r 0 v i ) _1 r 0 . (2.13) 

Now in Eq. (12.121) the factors before and after the square 
brackets refer to only one- and two-body interactions (the 
latter referring to interactions between bodies 1 and 2 , 
and 1 and 3, respectively), so the two-body interactions 
between 2 and 3, and three-body interactions are all con¬ 
tained in the quantity in square brackets. Now, in terms 
of the potential, the corresponding scattering matrix is 

T<=V i (l-r 0 V i )- 1 . (2.14) 

Introducing the modified scattering matrix defined by 
Shajesh and Scliaden [l3| . 

T = Tr 0 , (2.15) 

and using the cyclic property of the trace, we find the 
two and three body terms: 

£23 = ~ T¥ln(l - T 2 T 3 ), (2.16) 

which is sometimes called the TGTG formula, and 


E 123 — — — Trln ( 1 X 23 


X 21 T 2 (1+T 1 )X; 
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xf 3 (l + T 1 )-T 2 T 



(2.17) 


where 

X*j = (1 - Tif j) _1 . (2.18) 

Expression (12.171) has a rather evident geometrical in¬ 
terpretation in terms of multiple scattering off the three 
objects. This is not written in as symmetrical a form 
as in Ref. |l3l j , but is somewhat simpler, particular for 
the Casimir-Polder applications that follow, where body 
1 represents the atom, so is treated weakly. 


III. POLARIZABLE ATOMS BETWEEN 
PARALLEL CONDUCTING PLATES 

As a simple check of the machinery developed in the 
previous section, we revisit the calculation of the inter¬ 
action energy of an anisotropically polarizable atom be¬ 
tween parallel conducting plates, a geometry first ana¬ 
lyzed by Barton E2- Since we know the Green’s dyadic 
T for parallel plates, it is easy to derive the interaction 
energy from the general Casimir-Polder formula 

/ OO 

d( Tr a. • T, (3.1) 

-OO 

where the integration is over imaginary frequency, u> —> 
iC,. 1 Here the (in general, anisotropic) polarizability of 
the atom is a. In the following we will assume a is in¬ 
dependent of frequency, i.e., we are working in the static 
approximation. The interaction energy for one conduct¬ 
ing plate at z = 0, one at z = a, and the atom at z = Z, 
0 < Z < a, is 

an + a 22 — a 33 
Ecp = -W-^ 4 ) 

-^[C(4,Z/a) + C(4,l-Z/a)],(3.2) 

in terms of the Hurwitz zeta function. Here the two- 
body interactions between the atom and one or the other 
plate are isolated by extracting the parts singular as Z 
0 or Z —> a: 

C(4 ,Z/a) = ) + C(4,1 + Z/a), 

C(4,l-Z/o) = (^Tz) +C(4,2 -Z/a). (3.3) 


The total Casimir-Polder energy is the sum of two-body 
and three-body terms, 


Ecp = E \2 + E /13 + Ed 23 , (3.4) 


where 1 denotes the atom, 2 the lower plate, and 3 the 
upper plate. Here 


E12 — — 


tr a 
8 ttZ 4 ’ 


E i 3 = 


tr a 

87 t (a — Z ) 4 ’ 


(3.5a) 


1 This replacement requires knowledge of the analytic properties 
of the integrand. There can be serious subtleties involved in this 
“Euclidean transformation,” for example, see Ref. [24ll . See also 
Ref. [H|. 

2 The particular combinations of Hurwitz zeta functions occurring 
here and in the following are striking. Such combinations occur in 
several places, for instance when considering the Casimir energy 
for two parallel plates in D dimensional spacetimes, where the 
argument 4 is replaced by D. Two of the first papers in this 
direction are Refs. [26l .[27| . See also Ref. [2f|. Related structures 
appear for the theory of the piecewise uniform string [29l - [3l}l . 
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FIG. 1: Three-body contributions to the Casimir-Polder in¬ 
teraction of an anisotropically polarizable atom between two 
parallel polarizable plates. What is plotted is the ratio 
of the three-body contribution relative to the total energy, 
r = E123 /Ecp- The three-body contributions are gener¬ 
ally very small. They become appreciable only far from both 
plates, where the Casimir-Polder energy is very small. Plotted 
is the ratio for isotropic atoms (middle curve), atoms polar¬ 
izable only in the direction perpendicular to the plates (top 
curve), and polarizable only parallel to the plates (bottom 
curve). In the case when only a zz ^ 0, the sign of the three- 
body correction is the same as that of the two-body term. 


and 


El23 — 


ctn + a-22 — ol 33 


C(4) 


47TC1 4 

[C(4j 1 + Z/a) + £(4, 2 — Z/a)} . 

(3.5b) 


Note that the first term in £’123 is independent of Z, so 
it does not contribute to the Casimir-Polder force on the 
atom, but is a Casimir-Polder correction to the Casimir 
force between the plates. The two-body energies over¬ 
whelmingly dominate the Casimir-Polder interaction, as 
shown in Fig. [D For isotropic atoms, the largest three 
body correction is only a 0 . 8 % reduction at the mid¬ 
point between the plates, where the energy is very small. 
For atoms only polarizable parallel to the z direction 
the three-body correction is an 8 % increase at the place 
where the energy is the smallest, while for purely trans¬ 
versely polarizable atoms, the three-body correction is a 
6 % reduction at the midpoint. 


where R = (0, 0, Z) is the position of the atom. The free 
electromagnetic Green’s dyadic can be written as 

ro(i-r') = J ^e ik -(--')u 7o(Z)Z0) (3.7) 

where 


7 o(*,<) = (E + H)-U-'^-*'I, (3.8) 


with the usual abbreviation n = \/k 2 + 2 . Here E and 
H are matrices corresponding to the transverse electric 
(TE) and transverse magnetic (TM) modes, 

/ s 2 —cs 0 \ 

E = — ( 2 I — cs c 2 0 I , 

\ 0 0 0 / 

/ c 2 d z d z > csd z d z ' ikcd z \ 

H (z, z') = I csd z d z / s 2 d z d z ' iksd z I . (3.9) 

y — ikcd z ' —iksd z ' k 2 / 

Here k 2 = ky and c (s) is the cosine (sine) of the angle 
between the direction of kj_ and the £-axis, c = k x /k, 
s = ky/k. The polarization operators are transverse, in 
the sense that 


ik± • H + d z z • H = 0, (3.10) 

and similarly for E. Thus the modified scattering matrix 
for the atom is 

T 1 (r, r') = Ancx5{z - Z)5(r ± ) J 

x(E + H Mz,z')—e- K \ z - z '\. (3.11) 

2k 

The following composition properties of the E and H 
operators are easily checked: 

EH = 0, (3.12a) 

EE = -C 2 E, (3.12b) 

U(z,z')n{z",z'") = (, k 2 + d Z 'd Z ")H(z,z 

(3.12c) 

For a single plate, say a conducting plate 2 at z = 0, 
we have the reduced Green’s dyadic in the form 


A. Multiple-scattering calculation 


7 = E.g B + H 5 ff , (3.13) 


where 


Since in more general situations we do not have an 
exact solution available, we want to calculate the three- 
body corrections (I3.5bl) , using the multiple-scattering for¬ 
mula (12.171) . For this purpose, we need to compute the 
scattering operators for the three bodies. 

The scattering matrix for the atom is simply 


g E,H 


(z,z') = g 0 (z,z')T^-e- K (M+W 


1, 

sgn(~)sgn(2'), 


9 o(z,z') = (3.14) 

Then the reduced modified scattering matrix is 


Ti(r, r') = Vi(r, r') = 4naS(r — R)d(r — r'), (3.6) 


t 2 (z,z') = 7 0 1 (7~7o )(z,z'). 


(3.15) 
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This is evaluated by using the transverse property of E 
and H [hence, the vector Helmholtz operator reduces to 
—V 2 + C 2 —see Eq. and 



2 n8(z), 


(3.16a) 



sgn(z)e -K ^ = 25'(z)e- K W. (3.16b) 


Thus the modified scattering matrix for a conducting 
plate at z = 0 is 


*2 (z,z r ) = 


1 

C 2 L 


E (z,z')S(z)e-^ z '\ 


+ H(z, z')-5'(z)e-^sgn{z')e- K \ z '\ 

K 


(3.17) 


However, because 8 '{z) is an instruction to integrate by 
parts and evaluate at z = 0, which action is on the expo¬ 
nential propagators occurring in every case, we can use, 
as in Ref. [22l |. 

ta (z,z') = ^(E-H ){z,z')8{z)e-^‘\. (3.18) 


Then it is easy to see that the two-body interaction 
energy between the atom and the plate is as expected: 

E 12 = -TrTiTa = —(3.19) 
12 2 1 2 8 ttZ 4 V ' 


B. CP interaction between atom and two parallel 
conducting plates 

The three-body interaction is worked out by simpli¬ 
fying the multiple-scattering formula (12.171) for the case 
when there is only one interaction with the atom, since 
that coupling is always weak: 

E 123 = - TrX 2 3(T 2 T 1 T2T3 + T , 2T , 3T 1 T 3 

+ f 2 T 1 T 3 + f 2 f 3 T 1 ), (3.20) 

where the T operators are given in Sec. IIII Al 

Let us look at the E and H parts separately. For the 
TE part, 


X 23E = E T—s^Siz), (3.21) 


where integrating over the directions of kj_ gives for the 
trace 


tr a 



—> ( a n + a 2 2 )- 


(3.23) 


Thus the TE contribution is 

i /*C 

Oi 11 + 0-22 / 


pTE _ 
-^123 — — " 


127T 
-2 - 


dn n 6 


- 1 


g—2 kZ _j_ g—2«i(a— Z) 


an + a 22 


32?ra 4 
+ C(4,2 — Z/a) 


- 2C(4) + C(4,1 + Z/a) 


(3.24) 


The TM contribution is similarly worked out, with the 
result 


E- 


TM 

1 2 3 


1 f°° dKK 3 f an + a 22 


27T 


0 2 Ka 


- 1 


£ ~2kZ £ 2K(a-Z) 


2 

-2 


2 

+ g«33 


-2 kZ _j_ ^— c 1k{cl—Z') _|_ 2 


3(an + a 22 ) 

327ra 4 

x [—2((4) + C(4,1 + Z/a) + C(4,2 - Z/a)} 

- ^ [2C(4) + C(4,1 + Z/a) + C(4, 2 - Z/a)) . 

(3.25) 


Adding this to the TE contribution (13.24[) . gives the 
three-body energy (I3.5bl) . 

The three-body corrections are dominated by the 
three- and four-scattering contributions, given by the ex¬ 
plicit scattering terms in Eq. (13.201) . with the multiple- 
reflection quantity X 23 set equal to 1. That translates 
into replacing the zeta functions in Eq. (13.5bl) by their 
leading terms, £(4) —> 1, £(4, a;) —> 1/a; 4 . Figured] 
compares the exact three-body corrections to the leading 
three- and four-scattering approximations. (It is geomet¬ 
rically obvious why the odd-scattering terms give con¬ 
tributions which are independent of the position of the 
atom, because the path length is then an integer multiple 
of the plate separation.) 


IV. CASIMIR-POLDER INTERACTION 
BETWEEN ATOM AND WEDGE 


so 


pTE _ 
^123 — 


^2 Jd( (dkx)tra^-^ 


2kcl 


47 t 2 ./ V C 2 / 1 — e 2fta V 2 k , 

x e~ 2K,z + e - 2K(a " z) - 2 , (3.22) 


Our goal had been to compute the three-body correc¬ 
tions for an atom near an aperture created by two facing 
wedges, as shown in Fig. [3j Here we have two parallel 
conducting wedges, with opening angles /?, whose apexes 
are separated by a distance 2A'. As both interior wedge 
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z/a 


write this most conveniently in terms of the quantity 
k = \Jk 2 + ( 2 , where ( = —ico is the imaginary frequency 
and the wavenumber in the longitudinal direction or z di¬ 
rection, the direction perpendicular to the plane of the 
figure, is k. With 9 defined as the angle from the top sur¬ 
face of the wedge, that Green’s dyadic is (the prime on 
the summation sign is an instruction to count the m = 0 
term with half weight) 



FIG. 2: Three-body contribution to the Casimir-Polder in¬ 
teraction of an isotropically polarizabable atom between two 
parallel polarizable plates. The interaction energy is given in 
units of a/o 4 . The upper dotted horizontal line is the three- 
scattering approximation, the lower dashed curve is the four- 
scattering approximation, the upper solid (red) curve is the 
sum of these two contributions, which is only slightly above 
the the full three-body energy (solid blue curve). 




FIG. 3: Two facing conducting wedges, with an anisotropi- 
cally polarizable atom passing along a line perpendicular to 
the symmetry plane of the wedges. 


E(r, r') cos mpO cos mp9' 

+ H(r, r') sin mpO sin mp9' 
Imp (p'P< ) K fn p ( ) 


D ik(z — z f ) 


(4.1) 


Here p> (p<) is the greater (lesser) of p, p\ and p = 
— /3). The electric and magnetic polarization 
dyadic operators are 


E(r,r') 


H(r, r') 


- (p- p d e - 0d p ^j (p'^d e > - O'dp'J 

x(vi -k 2 ) 

-C 2 (z x Vj_)(z x V' x ) 

—V 2 (z x Vj_)(z x V x ), (4.2a) 

ik ^ pd p + 9-dg^j — zVy 


-ik ( p'dp' + 9'^-jC 


— zV 


(?fcV X — ZK 2 ) (-iW'l — ZK 2 ) 
[V x (V x z)][V' X (V' x z)]. 


/2 

X 


(4.2b) 


angles go to zero, the situation reduces to two conduct¬ 
ing half-planes, lying in the same plane, with a gap be¬ 
tween them. The Casimir-Polder interaction between the 
two wedges and an anisotropically polarizable molecule, 
located at coordinates p, 4> relative to the apex of one 
wedge, is to be computed. In particular, we wish to study 
the three-body interaction, which involves scattering off 
all three objects, as a correction to the more elementary 
calculation of the interaction between the atom and a 
single wedge, which is given in Ref. [4j. For the latter, 
repulsion can be achieved for the angle (f> sufficiently close 
to 7r, provided f3 is smaller than 108° = 1.88 radians. 


[Here the polarization operators differ from those given 
previously in Eq (13.91) by the replacements E —> —V X E, 
H —» — V X H. For a further discussion of the properties 
of the polarization operators, see Ref. [23).] In the sec¬ 
ond forms in Eq. (14.21) we have used the modified Bessel 
equation, that is, that for either modified Bessel function 


(-Vi + AC 2 )e ily((9 - e,) 


Iu{np) _ n 
K v (np) 


(4.3) 


In the following we will need the composition proper¬ 
ties of these operators, analogous to those in Eq. (13.121) : 


A. Scattering Matrix for a Perfectly Conducting 
Wedge 

In Ref. [|[ we gave the Green’s dyadic for a single 
perfectly conducting wedge, in terms of polar coordi¬ 
nates based at a point on the apex of the wedge. We 


E(r, r / )H(r // , r'") = H(r,r')E(r", r"') = 0. (4.4a) 
E(r, r , )E(r", r'") = E(r, - k 2 ) 

-» K 2 C 2 E(r,r'"), (4.4b) 

H(r, r , )H(r // , r"') = H(r,r- fc 2 ) 

^ K 2 C 2 H(r,r"'), (4.4c) 
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where we will understand that after differentiation, the 
intermediate coordinates r' and r" are set equal and in¬ 
tegrated over; that is, a spatial matrix multiplication is 
understood. 

To construct the modified scattering matrices, we need 
the free Green’s dyadic, which we can write as Eq. m , 
where a representation for the scalar Helmholtz Green’s 
function in cylindrical coordinates is 


Go(r-r') = 


dk ik(z-z') 1 

2tt 2tt 


OO 

E ■ 

m =—oo 


,1171(0 — 0') 


X Im (Kp<)A' m (Kp>). 


(4.5) 


It is easy to verify that in terms of the mode operators, 

r 0 (r, r') = — E(E + H)(r, r')G 0 (r - r'). (4.6) 

v x 


Using the above, we find the modified scattering matrix 
for the atom to be 

/ rlk 

— e M z ~ z ) 

„ 1 im(e-e'}Irn{KP<)K m {Kp>) 

2 tt ^ e -k 2 

m=—oo 

(4.7) 

To work out the T matrix for the wedge, we start from 
Eq. (12.91) . 


T = To Hr-r 0 ). (4.8) 

The inverse free Green’s function is the differential oper¬ 
ator given in Eq. (12.51) . It is easy to check that 

K 2 - V 2 

Eo X E = E 2 -H (4.9a) 

or 

K 2 - V 2 

r^H = h -(4.9b) 

CO 2 

The Helmholtz operator appearing as the last factor here 
would annihilate the scalar Green’s functions appearing 
in Eq. m, except on the boundaries, where the nor¬ 
mal derivatives give contributions to the scattering ma¬ 
trix that live entirely on the surface of the wedge. This 
is precisely the same as what occurred for the planes in 
Sec. 1III Al see Eq. (13.161) . Here, because we are consid¬ 
ering the region exterior to the wedge, we interpret the 
angular mode functions as follows: (9 £ [0, H]) 

cos mp6 —> cos mp9 r)(9) rj(fl — 9), (4.10a) 

sin mp9 —> sin mp9rj(9) ??(fl — 9), (4.10b) 

where H = 27 t — p is the exterior wedge angle, and the 
step function is defined by 


(4.11) 


Then we see that 


[dg + ( mp ) 2 ] cos mp9 = 5'(9) — (—1 ) m 5'(9 — H), 

[<9g + (mp) 2 ] sinmp9 = mp[6(9) — (— \) m 5(9 — fi)]. 

(4.12) 


From this we can immediately read off the modified scat¬ 
tering matrix for the wedge: 


Twedge (**, I” ) 



dk 1 
2-7T C / 2 


x jE(r, r') [S' (<t> - p/2) - (-1 rS'M + p/2)] 
x cos mp((j)' — p/2) 

+ mpH(r, r') - p/2) - (-1 ) m 6(</> + p/2)] 

x sin mp(<f> — P/2)Xe lk ( z ~ z ) 


X .2 2 Xnp(^P<)Rmp(^P>), 

K,p 


(4.13) 


where we have shifted to the angular variable </> mea¬ 
sured from the symmetry plane of the wedge, as shown 
in Fig.[3j p, (j)' £ [P/2, 2 tt — P/2], and the delta functions 
are understood to be periodically extended, with period 
27T. 


B. Two-body Calculation 


We now want to use this multiple scattering formalism, 
particularly Eq. (12.161) . to reproduce the results found in 
Ref. Q. Putting together the scattering matrix for the 
atom, Eq. (1X711 . and that for the wedge, Eq. (14.131) . we 
obtain the following expression for the Casimir-Polder 
energy, 


Enw — 


• P 4 . 

——r tr a. 
4tt 3 


d(dk f°° df/_ ^ y, 

K 2 Jo P' ^ , 

m=— oo m' = — oo 


x [mE(r, r") - m'pU(r, r")] e irn ' p ^"- 0/2) 


xe 


imcj) 


0 —im/3/2 


- (-l) m e 


m' impj 2 


X ^m {^P<)I^m (^P> )I\m' \p{^P< )-^|m / |p(^P>)’ 

(4.14) 


Here p<,> is the lesser, greater of p,p', and /5< i > is the 
lesser, greater of p ", p ’. After the differentiations con¬ 
tained in E and H are performed, the coordinates p" 
and p" are to be set equal to p and cf), respectively. 

Although we can carry out the m summation, or the p' 
integration, it seems difficult to bring Eq. (14.141) into the 
closed form given in Ref. Q. So we initially will content 
ourselves with a special case, P = tt or p = 1, that is, the 
interaction of an atom with an infinite conducting plane. 
In that case, we may as well ultimately set (/> = </)" = n. 
Then only Bessel functions of integer order occur, and 


v( x ) = 


1, x > 0, 
0, x < 0. 



















both the vn and ml sums can be carried out, using the 
addition theorem, 


OO 

Kq(kP) = (4.15) 

m =—oo 

where 

P = \Jp 1 + p' 2 — 2 pp' cos(</> — i ft). (4-16) 


Then the energy can be written as 
■ f°° dn 


Eap ~ 2tt 2 


f 0 


K 


/0 


^-trc*(E- H)(r,r") 


K 0 (k\/ p 2 + p' 2 - 2 pp' sin^>) 
1 d 


x --^ijj k o( k V p" 2 + P' 2 - 2 p"p' sin ft') 
i oft' 


1 

27T 2 


— (sin (j) —» — sin ft) 
/*00 

dn 


p"—>p, (f)"—>(j)=7T 

'OO 

dptr a(E — H)(r, r") 
p" cos ft' 


. yV 2 + p" 2 — 2p'p" sin ft' 
xKq(k(\Zp' 2 + p" 2 — 2p'p" sin. ft') 
xK 0 (Ky/p 2 + p' 2 — 2 pp' sin</>) 


+ (sin </>—}■ — sin 0) 


■ (4.17) 


p" p ,<t>" —*<{>= 7T 


The simplest situation occurs when the atom is only 
polarizable along the axis of the wedge, a = zza zz . Then 
tr aE = 0, tr aH = K 4 a zz , and we have 


F, = 

- LJ ap — 


_ r oo 

Otzz I j I (■ 

“ ./» P 7F 


p' 2 Jo 


dn. k 4 


xK 0 (k\/ p 2 + p' 2 )K\{k\] p 2 + p' 2 ) 


8irp 4 


(4.18) 


which is the expected Casimir-Polder energy. With only 
a bit more effort, we find the familiar result for arbitrary 
polarization, 


£cp = -0. (4.19) 

Given the difficulty of even analytically reproducing 
the two-body correction given in Ref. [dj, it is not sur¬ 
prising that we did not get very far with the three-body 
calculation. In the Appendix we consider the scalar ana¬ 
log for the two-body effect, and although we get a bit 
further, we have been unable to reproduce the analytic 
result obtained by the direct calculation. 3 So we turn, 



FIG. 4: An anisotropically polarizable atom, denoted by 
a, is symmetrically located relative to two identical parallel 
perfectly conducting cylinders (with axes coming out of the 
page). The centers of the cylinders are separated by a distance 
h, and each has radius a. The atom is on the line bisecting 
the line connecting the centers of the two cylinders, but a dis¬ 
tance y above it. The atom is a distance R from the center 
of either cylinder. The angle 6 is defined by cos 6 = y/R. 


instead, to another problem, the interaction between an 
atom and a pair of cylinders. 


V. CASIMIR-POLDER INTERACTION OF 
ATOM WITH TWO PARALLEL CYLINDERS 


The difficulties in extracting usable expressions for 
three-body effects for the atom-wedge-wedge problem has 
to do with the lack of a scale, so multipole expansions, 
for example, are not applicable. Therefore, we turn to 
another example, that of an atom interacting with a pair 
of parallel cylinders, illustrated in Fig. [4] Here it is as¬ 
sumed the two cylinders are identical, with radius a, and 
their centers are separated by a distance h. The atom 
is located on the line bisecting the line connecting the 
centers of the cylinders, a distance R from the center of 
each, and a height y above the centerline. We will con¬ 
sider the case when the atom is only polarizable along 
the bisecting line. 


A. Scattering matrix for cylinder 

We first need to compute the scattering matrix for the 
cylinder. The Green’s dyadic, given in Ref. ;32|], part 
of which appears in Ref. [5|, can be written, with slight 
notational changes, as 4 


r cy i(r,r') = - J2 


dk 


e im(4>-4>') e ik(z-z') 


m— — oo 


(2n) 2 

x \EF m (p,f/)+HG m (p,p')], (5.1) 


3 However, it is possible to recast the T 1 T 2 expression into the 
form of the direct CP energy Eqp = Tr f d£a(G — Go) —see the 
Appendix. 


4 As noted in Ref. [33| , only the terms in the Green’s functions 
involving modified Bessel functions, and not powers of the radial 
coordinates, contribute to the electric and magnetic fields. 
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where the TE and TM Green’s functions are, outside the 
cylinder, 


F'm (Pi P ) o 


Gm(p,p') = 


I'm ( K P < )K m (Kp > ) 


r m {na) 


K m (np)K m (Kp') 


K' m (na) 
dm («P<)-Km(«P>) 


Im(Ka) 

K m (na) 


K m (np)K m (Kp') 


,(5.2a) 


,(5.2b) 


and the polarization operators are the same as given 
in Eqs. (TOl) . The modified scattering matrix is given 
by Eq. (14.81) . where, because of the tranversality of the 
polarization operators, the inverse free Green’s operator 
may be replaced by Tq 1 —> ^-(V 2 — ( 2 ) [cf. Eq. (14.91) ]. 
which vanishes everywhere but on the surface of the 
cylinder. Because we have a perfectly conducting body, 
F rn (p, p') = 0 if p' > a > p, and so for the TE function 
with p' > a 


{F m -F 0 m ){p,p') = --K m {np') 


d'm( Ka ) 


K' m {na) 


xK m (np)r](p - a) + Im(np)r)(a - p) 


. (5.3a) 


Similarly, the TM functions are 


1 


( G m -G° m )(p,p') = -- 1 K m {Kp') 


Im(Ka) 

_K m (na) 

. (5.3b) 


xK m (np)r)(p - a) + I m (Kp)Tj(a - p) 

Then a simple calculation leads to the scattering matrix 

OO 

Tcyi = E 




m= — oo 


J ( 2^) 2 

„ 1 1 d K m (Kp') 

. npdp P P a K' m (na) 


k 2 C 2 a 


- Hd(p - a) 


K m (Kp') 


K m {na) 


(5.4) 


To check its validity, we reproduce the two-body inter¬ 
action between one cylinder and the anisotropic atom, 
for which we obtain 


El o — 


1 


27T 


TrT 


atom J- cyl 


a 


d(dk 

(27r) 2 J K* 

m =—oo 


E * 


4m («a) 
K' m (na) 


E(r,r') 


Im(na) 


H(r,r') 


E m (tm) 

xK m (K.p)K m (np') 


(5.5) 


r=r'=R 


This is the general result, which may be derived by sim¬ 
pler means. In particular, for the situation envisaged in 
Fig.® where R = h /(2 sind), and the atom is only polar¬ 
izable along the y direction, we obtain the formulas given 
in Ref. [5j 5 in terms of the distance of closest approach 
Ro = h/2, 


pill _ ol sin 4 9 ^ f°° j I m (Ka sin 9/Rq) 

CP “ 47 t i? 4 dxx Km{K asme/R 0 ) 

x \m 2 K^ n {x) sin 2 0 + x 2 K'£(x) cos 2 0] , 

(5.6a) 

te _ a sin 4 6( y- f°° I'mi^a sin 9 / Ro) 

CP “ 4 - R o dXX K^a S me/R 0 ) 

x \m 2 K^ n {x) cos 2 6 + x 2 K' 2 z (x) sin 2 6] . 

(5.6b) 

These formulas show that repulsion indeed occurs along 
the bisector (y) direction provided the distance of closest 
approach Rq is larger than about 7 times the radius of 
the cylinder; in this case the m = 0 term dominates, 
and the TM contribution is much larger than the always 
attractive TE contribution, as we will see below. 


B. Three-body correction 

We now want to see if the above repulsive effect sur¬ 
vives when the effect of both bodies are included. By 
virtue of the symmetry seen in Fig.® the two-body forces 
in the y-direction are doubled. That is, if the atom is 
called body 1, and the cylinders are 2 and 3, respectively, 
the two-body terms are just 

£"2— body = Ei 2 + E13. ( 5 - 7 ) 

The three-body terms are computed from Eq. (13.201) . In 
view of the remarks above, because we are considering 
large distances between the cylinders, it should suffice to 
consider the three- and four-scattering terms, the higher 
terms being suppressed, 

E 3 —body ~ E123 + E132 + E1232 + E1323, 

E123 = — — TrTiT2T 3 , E1232 = — — Tr T1T2T3T2. 

(5.8) 

Further, we expect dominance by the m = 0 TM mode. 
However, we can effect considerable simplification before 
we make that last approximation. Indeed, the formula 


5 The signs of the energies given there, and in Fig. 3 of that refer¬ 
ence, should be reversed. All physical conclusions in that paper, 
however, are correct. 
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and the corresponding TE 3-scattering energy is 
d( dk 1 


pTE _ 
^123 — 


(27r) 2 k 1 


■ tr aE(r, f) 


o2-k il/ 


E / Z 


i(m<f>—m'4>) i(m'<fi'—m<f>') 


27 T 


x Fr7§Tir' f -<^-'M 

A m(« a ) K L'( Ka ) 


P=P—R 

(5.10) 


where a and </>' are the cylindrical coordinates of a point 
on the surface of the second cylinder relative to the cen¬ 
tral axis of that cylinder, the same point being located 
at cylindrical coordinates p' and <fi' relative to the central 
axis of the first cylinder. These coordinates are related 
by 


FIG. 5: The Casimir-Polder energy of an anisotropic atom 
or nanoparticle passing on the symmetry line perpendicular 
to a pair of identical perfectly conducting parallel cylinders. 
The energy, apart from a factor of a/(47ri?o)i is plotted versus 
<j> = 7 t /2 — 9, that is, the angle between the line connecting 
the axes of the cylinders and the line connecting the atom 
with the center of either cylinder. The bottom curve (red) is 
the two-body energy, the top curve (dotted-blue) is the three- 
scattering correction, the second-fronr-top curve (dashed red) 
is the four-scattering term, while the second curve from the 
bottom (blue) is the total Casimir-Polder energy. The ener¬ 
gies are plotted for a/Ro = 0.01. Because a/Ro 1 it is 
sufficient to include only m = 0 for the three-body correc¬ 
tions. Also, the TE contributions are completely negligible. 
It is seen that the three-body effects, in fact, are very small, 
and do not significantly alter the repulsion between the atom 
and the pair of cylinders, and the four-scattering terms are 
quite negligible. 


i r, . ,, a sin d>' 

p = h 2 +a 2 — 2ahcos(f> , tan^> =-=-. (5.11) 

a cos </>' — h 

The atom is located in the two coordinate systems at 

R = (Ri <£, 0 ) = (R, 4>, 0 ), (5.12) 

where 

^=| + 6 '> ( 5 - 13 ) 

It would now be straightforward to work out the multi¬ 
pole expansion of Eqs. (15.91) . (15.101) that is, a power series 
expansion in powers of a/h. We will content ourselves 
with the lowest term, which means we can set m = 0 , 
because only for small values of a/h do we have 2 -body 
repulsion. Higher terms are suppressed by powers of a/h. 
In this limit, the TE energy is completely negligible, be¬ 
cause of the appearance of derivatives of Bessel functions. 
The behavior of the Bessel functions for small argument 
makes this point clear: 

I 0 {z)~l + J z 2 , K 0 (z )-(5.14a) 

(-) ~ k o(z) - z->0. (5.14b) 

2 Z 


Then we get a very simple explicit formula, again for an 
atom polarizable only in the bisector ( y ) direction, 

E™ „ _/^yy 2 q in 4 q 

123 47rl?g 

x f°° dx x 3 K ° Vlr Sm 0)1,1 (xa si 111 0/Bo 1 K' 2 (x ). 
Jo Kfi^xa sin Q/Rq) 

(5.15) 

Here h = 2 Rq is the separation distance between the 
axes of the cylinders. In the same approximation, the 
4-scattering contribution is also given simply: 

pTM _ a VV 

jC /1 oqo i cos u sm u 

1232 47ri?o 

x r dx x 3 A ° 2(2x Sin g)/ ° {xa Sin d/R ° ] K' 2 (x). 
Jo K[/(xa sin 9/R 0 ) 

(5.16) 

These TM corrections are plotted in Fig. [5] as a func¬ 
tion of <j> = 7r/2 — 9, so 4> = 0 at the position of the 
atom closest to the cylinders, and compared to the two 
body contributions. (All the TE corrections are com¬ 
pletely negligible.) The figure shows the two-body CP 
energy, entirely dominated by the m = 0 TM contribu¬ 
tion, and the m = 0 three-body correction, dominated by 
the three-scattering terms, in the limit of large distance 
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between the cylinders. The TE contributions were com¬ 
puted, but found, as expected, to be completely insignif¬ 
icant in the large distance regime. The TM three-body 
correction is not quite negligible, but does not affect the 
Casimir-Polder repulsion discovered in Ref. [&]. 


VI. DIFFICULTY OF ACHIEVING HIGH 
ATOMIC ANISOTROPY 

In order for Casimir-Polder repulsion to be possible, 
the atom interacting with a body must have a sufficiently 
anisotropic polarizability tensor. Defining an anisotropy 
factor 7 according to 

a = a 22 zz + 7 « 22 xx + 7 a 22 yy, (6.1) 


implying 


7 = 


Tr(a) - a 22 
2 a 22 


( 6 . 2 ) 


we found in Ref. [4] that for an anisotropic atom inter¬ 
acting with a half-plane, the critical value of 7 was 1/4. 
Values 7 < 1/4 give repulsion in certain circumstances, 
whereas for 7 > 1/4 no repulsion is possible. The non- 
retarded interaction of an atom and a circular aperture 
was considered in Ref. [34j and the critical value of 7 was 
found to be 1/4 also in this case Q. 

In this section we investigate what minimal value of 
the anisotropy parameter can be achieved by preparing 
an atom in an excited eigenstate | nlm). Here, n denotes 
the principal quantum number, l = 0 ,...n — 1 is the 
quantum number for the orbital angular momentum and 
to = —/.../ that for its ^-component. The question is 
of great interest especially in the light of recent advances 
in experimental techniques using Rydberg atoms, atoms 
excited to high principal quantum numbers, near bound¬ 
aries (3f| HH, and noting that Rydberg atoms can take 
highly anisotropic shapes. 

Because of the close spacing of energy levels for highly 
excited states, and the fact that transitions to the states 
nearest in energy to | nlm) dominate the CP energy [36| . 
the Casimir-Polder interaction of a Rydberg molecule is 
essentially non-retarded even at atomically large separa¬ 
tions, up to hundreds of micrometers. It was shown that 
in such cases the interaction energy is proportional to the 
atomic dipole moment tensor [371 . 138l |, 


E nlm = -<dd> : VVG 


= -<dd> : r w= 


= 0 , 


(6.3) 


u—0 


according to Eq. & For convenience we will work with 
the ratio 


q = (*„)/( d 2 )- 

Defining 7 similarly as before 

(d 2 ) - (<4) _ 1 n 

m z ) 2 \q 



(6.4) 

(6.5) 


we will consider energy eigenstates | nlm) such that 
anisotropy becomes maximal, i.e., 7 becomes minimal 
and q maximal. 

To evaluate the anisotropy parameter, we insert the 
completeness relation \n'l'm'){n'l'm'\ = /, 

= E n'i'm' ( nlm\d z \n'l'm') {n'l'm'\d z \nlm) 

The dipole-matrix elements can conveniently be calcu¬ 
lated by means of the Wigner-Eckart theorem {39j, [4f| 


( n'l'm'\d s \nlm) = (— 1 ) 


l—m 


v 1 1 


— TO S TO 


(n 1 1'\\d\\nl) 


(6.7) 

where (nT||eZ||nZ) denotes the reduced matrix element 
and the Wigner 3 -j symbol can be given in terms of 
Clebsch-Gordan coefficients as 41] 


hi J2 j\ j | ■ _ m) (68) 

\mi to 2 mj ^2j + 1 

Substituting these relations into Eq. (16.61) and using the 
orthonormality relation [4lj 


V 2 3 + 1 X! 

777-17712 


( i 1 

32 

h 

( h 

32 

•M 

V TOi 

to 2 

m J 

V TOi 

to 2 

m J 


we find 

qim = \ 

l' 


&jj' ^mm/ 5 

(6.9) 

( 6 . 10 ) 


As expected from the symmetry of the problem, the 
anisotropy parameter depends neither on the reduced 
matrix element nor on the principal quantum numbers. 

The Clebsch-Gordan coefficients in Eq. (16.101) can be 
evaluated explicitly, leading to [tO] 


l 2 - TO 2 ( (/+ l ) 2 - TO 2 (ciM 

qim _ (2Z-1)(2Z + 1) + (21 + 1)(2Z + 3) ‘ ^ ‘ ’ 

For a given l , the anisotropy parameter obviously takes 
its maximum value for m = 0 , 

l 2 (1 + 1) 2 2Z(/ + 1)-1 

®° “ (2l-l)(2l + l) + (2l + l)(2l + 3) ~ 4Z(Z + l)-3 ' 

( 6 . 12 ) 

The latter expression is equal to 1/3 for / = 0, approaches 
1/2 for / -* 00 and takes it maximum value 

Qio = | => 7=3 (6-13) 

for 1 = 1. The maximally anisotropic eigenstate of orbital 
angular momentum is thus a p state. 

Since qi m is positive for any given choice of quantum 
numbers, it immediately follows that the anisotropy pa¬ 
rameter 7 is bounded below by 1/3 for any incoherent 
superposition of energy eigenstates. It is possible that 
stronger anisotropies could be realizable with a coher¬ 
ent superposition of states. However, the more likely 
venue for discovering such repulsive effects would be with 
anisotropic particles, such as elongated needles. 
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VII. CONCLUSIONS 

One of the principal features of Casimir or quantum 
vacuum forces is that they are not additive. Unlike clas¬ 
sical electrodynamics, one cannot simply sum pairwise 
forces. Such approximations clearly are invalid even for 
the simplest situations of parallel plates. This, of course, 
makes calculations more challenging. 

In this paper, we have explored some aspects of three- 
body interactions in the context of Casimir-Polder forces 
between an anisotropically polarizable atom or nanopar¬ 
ticle and two conducting surfaces. First we examined 
the role of such forces involving an atom between two 
perfectly conducting plates, a well-known problem 53, 
but one in which we can test our formalism and isolate 
explicitly the three-body terms. Then we turned to the 
interaction of such an atom with a pair of wedges; we 
reproduced the repulsive effects seen for an atom inter¬ 
acting with a single wedge j|J, but a closed form for the 
three-body correction remains elusive. So we then exam¬ 
ined the interaction of an anisotropic atom with a pair of 
parallel cylinders. The two-body repulsive effect found 
earlier [5] was reproduced, and the three-body correc¬ 
tion was computed in the limit of large separation be¬ 
tween the cylinders, which is the regime where repulsion 
is expected. The three-body correction is non-negligible 
in this limit, is completely captured by the TM tlrree- 
scattering approximation, but is too small to affect the 
earlier-found repulsion. 
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Appendix: Scalar analog for atom-wedge problem 


In this Appendix we will consider a scalar analog of 
the atom-wedge problem. Let the atom be described by 
the potential 


V/tom = 47 ra< 5 (r - R), (A.l) 


where R = (R , </, 0) is the position of the atom. The 
modified scattering matrix for the (Dirichlet) wedge is 
[note that the sign is reversed compared to Eq. (14.81) ] 


7 wedge — 1 Gq G w , 


(A.2) 


where 

G w { r,r') = BE f°° £ sinrn p(cf> - /9/2) 

x sinrop(<// - j3/2)I rnp (np < )K mp (np > ). 

(A.3) 

Applying Gg 1 = —V 2 + ( 2 we obtain the scattering ma¬ 
trix on the wedge: 


^ 7-oc P 

— (—l) m (5(</ + j3/2)\ smmp(q / — (3/ 2) 

(A.4) 

The two-body energy is given by 


P i fdoj ~~ i r ck 
E a = 2 / - 


Tr V[GqT 2 . 


(A.5) 

If we use the two-dimensional representation for the free 
propagator, 


r " 1 


2 „'.S A '" Wr±l) ’ 


Go(r) = 

the two-body energy can be written as 


(A. 6 ) 


pOO poo 7 / °° 

E 12 = ——^p 2 / dn.K / -j- ^ m sin mp{<j) — (3/2) 

n Jo Jo P m=1 

X Imp {^P< ) E mp (ft/?> ) 

x[K 0 (kP+) - (-l) m K 0 ( K P_)], (A.7) 

where P± is the distance between the atom and a point 
on the upper (lower) wedge boundary, 

P± = sjFt 2 + p' 2 - 2Rp' cos(</> =F (3/2). (A. 8 ) 

Now the integral of the three Bessel functions can be 
performed: 


dt t I v ((t) AY (t) K 0 (P± t/p>) 


1 ^ ^ +n+ 1 sin(re+ !)(</>-ff/ 2 ) 


2£ sin(</> — /3/2) 


v + n + 1 


(A.9) 


where ( = p < /p > . The radial integrals are then easy, 
and we are immediately led to 


E\2 = 


2 00 

■^ 2^2 m sill mp{(/) - 13/2) 

m =1 


y _-_ 

^ (n + mp)(n + mp + 2 ) 

n =0 

sin(n + l)(</> ~ 13/2,) 
sin (cj) — (3/2) 

_ , n m Sin(u + l)(/) + /3/2) 
sin((/ + (3/2) 


(A.10) 
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Now we replace the sum over m by an integral, 


E 


i 

mp + N 


gimpO 


00 /> 00 

53 y dt e -^ mp+N) e irnpe 

.(A.ll) 


m=l 

poo 


dt e 


—tN 


1 


e p(t-i8) _ l ' 


Then the n sum can be carried out as a geometric series, 
and the result is a single integral, 


-Ei 2 — — 


ap 


Att 1 2 R 2 


dt sinht sinhpt sin(</> — (3/2) 


1 


[coshpt — cos p(4> — /3/2)] 2 cosht — cos(</> — /3/2) 


_l___1_ 

[coshpt + cos p(<j> — (3/2)] 2 cosht — cos(</> + /3/2) J 

(A. 12) 


Alternatively, we can directly calculate the two-body 
energy from 


1 dC 

El 2=2 j ^^(G-Go). ( A - 


13) 


which may be directly evaluated in closed form: 

" 2 1 


Ei 2 = — 


8nR 2 


P 


svcf p{<j) — (3/2) 3 


+ o(W) 


. (A.14) 


The form (IA.12I) may be evaluated straightforwardly in 
analytic form for the two special cases p = 1 {(3 = n, that 
is, a plane), and p = 1/2 (/3 = 0, that is, a half-plane), 
in agreement with Eq. (1A. 141) . 


Et2{p=1) = - ^-L- , (A.15a) 

£ '^= 1 * / 2 > = ^fe + 1 ) <A ' 15b) 


For other values of p the analytic evaluation of Eq. (IA. 121) 
seems nontrivial; however, the integral is rapidly conver¬ 
gent, and the numerical coincidence with the closed form 
(1A. 141) is easily verified. 

Not surprisingly, it is possible to show that the explicit 
form of Eq. (IA. 131) follows from the multiple scattering 
form of the two-body energy (IA.7D , written in terms of a 
sum over four Bessel functions and a sum over m and m! 
as in Eq. (14.141) . This involves using the integral [notation 
as in Eq. (14.141) . with v = m'p], 


dp' 


-Im (np<) K m (kp > )I v (Ep< )K V (np> ) 
-\K v {kR)I v {kR) - I m (KR)K m (nR )\, 


(A.16) 

where the resulting two terms in the energy are summed 
over m and m ', respectively, using 


E 


gimQ 

m 2 — v 2 


m ——00 


777 / p sin m p 
^ m 2 — (m'p) 2 

m '=1 v ' 


E (- 1 ) 7 "' 

m '=1 


m'p sin m'p 
m 2 — (m'p) 2 


7 r cos v(6 — 7 r) 


v sin 7 tv 


(A.17a) 


7r . to mn 

-sin — (it — pO) esc-, 

2 p p p 

(A.17b) 


7r . „ TO7T . 

-smmtlcsc-. (A.17c) 
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